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Theory for Axial Rigidity of Structural Members 
Having Ovaloid or Square Perforations 

By Martin Greenspan 



Formulas are derived for computing the over-all lengthening (or shortening) of a tension 
(or compression) member having a uniform gross cross section and a series of similar perfora- 
tions of approximately ovaloid or approximately square shape uniformly distributed along 
the length. 



L Introduction 



In a previous paper [1] l the author demon- 
strated a method for computing the axial rigidity 
of a tension or compression member having a 
uniform gross cross section and a series of similar 
perforations of circular or elliptical shape uni- 
formly distributed along the length. The axial- 
rigidity factor K is defined so that KEA g is the 
rigidity that should be used in place of EA g in the 
ordinary formula for computation of the extension 
of the member. Here E is Young's modulus of 
elasticity, and A g is the gross cross-sectional area. 
The present paper extends this method to other 
shapes of perforation, originally described in 
reference [2]. The boundary of the hole has the 
parametric equation 

x=p cos fi-\-r cos 3/3, y=q sin fi— r sin 30. (1) 

Equation 1 represents a closed curve having 
symmetry about the #-axis and about the y-axis, 
and which, for appropriate values of p, q, and r, 
does not cross itself. By adjustment of the 
values of p ) g, and r, a variety of curves is obtained, 
including a good approximation to an ovaloid (a 
square with a semicircle erected on each of two 
opposite sides) and a good approximation to a 
square with rounded corners, as well as exact 



i Numbers in brackets indicate literature references at the end of this paper. 
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ellipses (r=0) of any eccentricity. The approxi- 
mate ovaloid obtained by taking 



p=2.0G3, g=1.108, r=-0.079, 



(2) 



is shown compared to the actual ovaloid in figure 
1. The approximate square obtained by taking 



p=2=l, r= 



-0.14 



(3) 



is shown in figure 2. The sides of the square 
are parallel to the axes of coordinates. By taking 



p=q=l, r=0.14, 



(4) 
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the same square, but with the diagonals parallel 
to the axes of coordinates, is obtained. The radius 
of curvature at the midpoint of the fillet is about 
0.086 times the length of the side of the square. 



It is shown in [1 , eq 2] 2 that for a member having 
perforations spaced 21 apart along the z-axis 
(parallel to the load), the axial rigidity factor K 
is given by 



1=1+4 



I u p dydz } (5) 




Figure 1. — Actual and approximate ovaloids. 

The dashed line represents the actual ovaloid and the full line the approx 
imate ovaloid of equations 1 and 2 



in which u p is the ^-component of the displace- 
ment of the boundary* of the hole, and the surface 
integral is extended over half the boundary of 
the perforation. The values of the load P and 
the displacement u p may be assumed, for suffi- 
ciently large values of I and of A g compared to the 
dimensions of the perforation, to be the same as 
in an infinite plate under uniform stress S x dis- 
turbed by a single perforation. 

For a member having circular holes, and of 
gross cross-sectional area A g and net cross-sectional 
area A nj 



P=S x A g C(n) =S;A(i-^-^)> 



(6) 
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where S x is the mean stress on the gross cross- 
sectional area and n=A g /(A g — A n ). Equation 6 
is derived in [1]. C(n) may be considered to be 
a correction factor arising from the finite area 
of the member. It has been shown [3] 3 that 
C(n) is nearly independent of eccentricity in 
case of an elliptical hole, and since it is nearly 
unity in any practical case, eq 6 may be used 
without serious error for a hole of any shape. 

The displacement u v may be calculated from 
the stress function given in [2] or by the method 
of Mushelisvili. 



Figure 2. — The approximate square of equations 1 and 8. 



* This notation is used throughout for the numbered equations in the 
references. 
1 See p. 545 of reference [3], 



II. Displacements 



1. Displacement-function method 

This method 4 presupposes a knowledge of the 
Airy stress-function, <t>. The form of is exhibited 
in [2] in terms of curvilinear coordinates (a, fi) 
such that eq 1 of the boundary of the hole reduces 
to a=a . The parameters a, t, and c of [2], 
which, together with a , define the shape of the 
hole, are related to the p, q, and r of eq 1 by 

e*=2±2, abe- a o= p -^, ac z e~^=r. (7) 



The function \p defined, in the absence of body 
forces, by 



6V 6 2 



dxdy dx 2 dy 



— 2~ V <P, 



(8) 



I See p. 130 of (Si 



1 3^ M - N d<t> 



(9) 
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with \[/ adjusted by means of the arbitrary func- 
tions of integration (of the form j t (x)-\-f 2 (y)) so 
that VV=0, is called the displacement function. 
The ^-component of the displacement is given by 



where /z is Poisson's ratio. In the present case the 
only values of u required are u P , those on the 
boundary of the perforation. On this (free) 
boundary, d<t>/dx is constant, and hence zero, since 
<j> is even in x and in y. Hence from eq 9, 



EUv=z lhj ( a==a °)' 



(10) 



According to [2, eq 18], 

<t> = C\ fa + C 2 fa + C%fa-\- Gtfa + 

0$4k + C 6 6 + Cr^7 + Os^s- 

Of these terms, all but those in fa, fa, and fa are 
harmonic and need not be considered for purposes 
of eq 8. Further, C 6 =0 for no shear, as shown by 
[2, eq 23]. Thus only the terms in fa and fa are 
required. From paragraph 5 of [2] fa and fa are, 
omitting harmonic terms, which disappear for 
operation V 2 , 

fa=4:ye a sin 0— 4afo/6~ a sin p—x 2 —y 2 , 

and 



fa=2ye~« sin 0, 



so that 



V^i^St-^ sin p—abe~ a sin 0)— 4, 



and 



V 2 2 =-4^-«sin 0. 



Equation 8 becomes 



= (7] %-(e a sin p—abe~ a sin /3) — 4 



a 

2 &y 



— 4C2*r-e~ a sin 0, 



from which 



^ = Ci 8(6 a sin p—abe-« sin 0)— 4y 

-4<7 2 e-« sin 0+Ifc+#, (11) 

in which Dx+E is the most general function of 
integration of the form j x (x) which preserves the 
harmonic character of d\p/dx. Because \f/ is har- 
monic by definition, d\f//dy as well as d\p/dx is 
harmonic, and application of the Cauchy- 
Riemann test to the expression 

dy dx 

shows it to be a function of the complex variable 
x+iy. Thus d^/dy is the real part of the function 



of which dypjdx is the imaginary part, so that from 
eq 11 



oy" 



--Ci(Se a cos $+&abe- a cos /3 — 4a;) 



+4<7 2 6-« cos p—Dy+E'. 

The constants D and E l are obtained by substitut- 
ing for d\f//dy in eq 9, which gives 

Eu=Ci(8e a cos p+8abe- a cos 0— 4x) 

+4C 2 e-° cos p-Dy+E'-(l + v)^ x - 

For ar=0(j8=±W2), ^=0 for all y, and d<j>jdx = 
by symmetry. Therefore, D=E' = Q, and 



6ty_ 

"9y " 



=Ci(86 a cos p+ Saber a cos p—4x) 



+4C 2 e-« cos 0. (12) 

Substitution of eq 12, with C\ and C 2 replaced 
by their values from [2, eq 23] into eq 10 gives 



Eu p 
& 



■■—x - 



2e aQ cos g 
'l-ac 3 <r 4a ° 



(2-<zcV 4a °-a 2 &c 3 <r fiaf °) ? 



where # is the value of x on the boundary a = a , 
which, together with eq 7, gives 



Eu 



(p-\-q) 2 —2pr 



£_ ~ _i_9 

S x °^ p+q-2r 



cos 0. 



(13) 



2. Mushelisvili's method 



This method, which is treated in detail by 
Sokolnikoff [4], expresses the stresses, displace- 
ments, and the boundary conditions in terms of 
two analytic functions, <£ and \p (not the same as 
those of the preceding section), of the complex 
variable z=x-\-iy. The functions and \p, which 
satisfy the problem of [2], have been given by 
Morkovin [5], who uses the notation 



and 



«=(P+ff)/2, *=(p-2)/2, r=r, 



so that the function 



(14) 
(15) 

(16) 



maps the boundary a=a onto the unit circle 
|f| = l. The value of f on the unit circle is 
denoted by a; thus o 7 , the conjugate of <r, is 1/a. 
The notation 

^(2)=*[«(r)]=*(rt, 

is used, and similarly for yf/. 
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The functions </> and \p as found by Morkovin 



are 



*(r)=«Br+£--a 






+s^+iT xy )^~a l - s t+Br^+Bi^l+~^ f 



where 



a v =Bt 



T + 8 , B's 2 , iT, ?/ S 2 



and 



r— s r—s r+s 

4B=S x -{-Sy, 2B'=S y —S x . 
In the present case, S v = 7^=0, giving 

*0)-f(tf+f-j*} (17) 

^)-&[-(i +(f+ ,,)($^±|) 

(18) 



The rectangular components, u and ^>,fof the dis- 
placement are expressed in terms of the functions 
(j> and yp by means of the relation 

E(u+iv) = (3 -,)*«•) - (1 +,) [^^+7(1)1 i 
so that on the boundary f =a, 

-(i+^r^+^i. (19) 

Substitution of the functions co, <j>, and ^ from 
eq 16, 17, and 18 in eq 19 yields, after considerable 
reduction, 

E fiA _.\ , , 2s 2 -2rt 1 

where 2 is on a = a , and therefore, 

fia, , 4s 2 -2sr-2rt _ '* ^ 

-= — ZoH cos /3, (20) 



& 



s—r 



-|-2.sr-yf+rr 3 + 



where 



a = 



4aj _ J (r+s)— 2s 2 
£* ~ r—s 



since o-=£ i/3 from eq 15. 

In terms of the original constants (eq 14), eq 20 
becomes, 

which is the same as eq 13. 



III. Axial Rigidity 



1. General 



Substitution of P from eq 6 and u p from eq 13 
into eq 5 gives 

'k= 1+ jghd.~J xdy 

+ 2 ^^ r J cos H' (21) 

where t, the thickness of the plate in which the per- 
foration occurs, is fdz. The two line integrals in 
eq 21, which are to be taken along half the bound- 
ary a=a , are obtained from eq 1, as 

J x(hj^A u -=^ (pg-3r 2 ), J cos pdy=^, 

where Aq is the area of the perforation. 

Equation 21 becomes 

</(P + 2)(P + 2<z) 
l-i i V* -2 M r+3r 2 (p+q-2r) 
K~~ L ^V g C(n) (p+%-2r)(pq-3r 2 ) ' {ZZ) 
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where V =A t is the volume of the perforation 
and V g =2A g l is the gross volume of one bay of 
the member. 

For an ellipse, p=a, q = l, r=0, where a and b 
are the semiaxes (a parallel to the load) and eq 
22 reduces to 

*r- 1+ coo tv 

which is the same as [1, eq 48]. 

It remains to compute the coefficient of 
Vo/V g in eq 22 for the various cases of interest. 

2. Ovaloid; load parallel to long axis 

Here 2>=2.063, £=1.108, r= — 0.079 (eq 2). 
Equation 22 becomes 



1 2.048 Vp 

K-^C(n)V g ' 



(23) 



In practical numerical cases, eq 23 gives sub- 
stantially the same results as [1, eq 13], derived 
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by a more approximate method. However [1, eq 
13] may be used for an ovaloid of which the rec- 
tangular base is not square whereas eq 23 applies 
only to the shape shown in figure 1. 

3. Ovaloid; load parallel to short axis 

In this case, p=1.108, (7=2.063, r=0.079. 
Equal 23 becomes 

* !, 4-968 Vo 

This case and the following cases are not covered 
by previous work. 



4. Square with Rounded Corners 

For the approximate square of eq 3 and figure 
2, p=q=l, r^— 0.14. Eq. 22 becomes 



*->- 



2.989 V 
' C(n) V e ' 



(25) 



for load parallel to the side of the square. 

For load parallel to the diagonal of the square 
(eq 4) p=g=l, r=0.14, and eq 22 becomes 



*-'- 



3.596 Vo 

" c(n) v; 



(26) 



IV. Summary 



Approximate formulas have been developed for 
the computation of the axial rigidity of a long 
tension or compression member containing a 
plate of constant thickness uniformly perforated 
with a series of similar holes of various shapes. 
The axial rigidity factor K is defined so that 
KEA g is the rigidity that should be used in place 
of EA g in the ordinary formula for computation 
of the extension of the member. 

In many cases of interest the axial rigidity 
factor K is given by an equation of the form 



1 _, , / Vo 



(27) 



where 

n=A g /(A g —A n ), A g being the gross and A„ 
the net cross-sectional area of the member 



ow-i-£.- 



1 

2n 4 



V =ihe volume of the perforation 
V g =tfo.e gross volume of one bay of the member 
jf=a constant depending on the shape of the 
perforation and the direction of the 
applied load. 

Values of the constant / for various cases are 
given in table 1. 

Table 1. — Values of f in equation 27 



Perforation 


Load parallel to- 


f 


Ellipse of semiaxes a and b . . 


Major axis, a 


l+(26/a) 


Do 


Minor axis, 6 


l+(2o/ft) 


Ovaloid (fig. 1) 


Long axis 


2.048 


Do--_ 


Short axis 


4.968 


Square (fig. 2) 


Side 


2.989 


Do _... 


Diagonal... 


3.596 
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